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ABSTRACT 


Imaging  spectrometry  requires  precise  knowledge  of 
wavelength  in  order  to  perform  various  spectral  analyses. 
This  thesis  project  was  tasked  with  performing  wavelength 
calibration  as  part  of  the  ongoing  development  of  the  Naval 
Postgraduate  School's  Lineate  Imaging  Near  Ultraviolet 
Spectrometer  (LINUS) .  This  calibration  was  necessary  for 
the  ability  of  LINUS  to  detect,  to  classify,  and  to 
quantify  several  different  chemical  species  by  observation 
over  the  near-ultraviolet  wavelength  band  of  200  to  400  nm. 

Experiments  were  conducted  to  detect  diffracted 
emission  lines  from  a  platinum  hollow  cathode  lamp  by  using 
the  LINUS  optical  train  and  five  different  UV  filters.  A 
Matlab  program  was  developed  to  compare  the  catalogued 
wavelengths  of  this  known  ultraviolet  source  with  the 
emission  line  positions  observed  on  the  LINUS  detector. 
Affine  transformation  and  cross-correlation  of  the  data 
produced  wavelength  calibration  curves  for  the  LINUS 
detector  in  each  of  the  five  associated  UV  wavelength 
regions . 
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I .  INTRODUCTION 


A.  PROJECT  CONTEXT 

The  completion  of  this  thesis  project  marked  the 
passing  of  another  milestone  in  the  continued  development 
of  the  Naval  Postgraduate  School's  Lineate  Imaging  Near 
Ultraviolet  Spectrometer  (LINUS) .  Upon  becoming  fully 
operational,  the  LINUS  instrument  will  be  the  third  imaging 
spectrometer  designed  and  constructed  at  NPS  and  will 
possess  greatly  improved  capabilities  over  its  forerunners. 

The  study  of  remote  sensing  has  produced  significant 
experience  with  imaging  spectrometers  in  the  NPS  Physics 
Department.  The  first  such  instrument,  the  Dual-Use  Ultra- 
Violet  Imaging  Spectrometer  (DUUVIS) ,  was  completed  in 
1996.  In  1997  a  follow-on  instrument,  the  Naval 

Postgraduate  School  Ultra-Violet  Imaging  Spectrometer 
(NUVIS) ,  possessed  notably  improved  capabilities  over  its 
predecessor.  NUVIS  demonstrated  the  ability  to  accurately 
detect  and  quantify  SO2  emissions  present  in  distant 
chemical  plumes  (Marino, 75).  Following  the  successful  use 
of  NUVIS  to  detect  a  single  potentially  hazardous  chemical 
species,  the  LINUS  project  was  conceived.  LINUS  was 
designed  for  significantly  enhanced  capability  over  NUVIS, 
especially  with  regard  to  the  measurement  of  spectral 
images  over  a  much  larger  wavelength  band  at  improved 
sensitivities. 

B .  PROJECT  OBJECTIVE 

Spectral  image  data  from  the  LINUS  instrument 


of  intensity  measurements 


as  a  function  of 
1 


consists 

spectral 


wavelength  and  of  two-dimensional  spatial  coordinates 
within  an  image.  In  order  for  wavelength  data  to  be  useful 
for  analyses  such  as  chemical  identification  via  spectral 
absorption,  an  accurate  wavelength  calibration  of  the 
instrument  must  first  be  performed.  Furthermore,  the 
capability  of  the  LINUS  instrument  to  observe  data  over 
several,  perhaps  discontiguous  subdivided  spectral  regions 
necessitates  the  requirement  for  a  rapid,  pseudo-autonomous 
method  for  wavelength  calibration.  Therefore,  this  thesis 
project  was  developed  with  the  goal  of  developing  a 
software-based  method  for  wavelength  calibration  utilizing 
known  wavelength  standards  and  the  LINUS  instrument  in  the 
non-imaging  spectrometer  mode  of  operation. 


C .  OUTLINE 

This  thesis  is  organized  into  six  chapters.  After 
this  brief  introduction.  Chapter  II  discusses  the 
fundamentals  of  hyperspectral  imaging  and  gives  an  overview 
of  the  LINUS  design.  Chapter  III  outlines  the  wavelength 
calibration  process  and  details  the  experimental  setup. 
Chapter  IV  describes  the  mathematical  calculations 
necessary  for  automated  wavelength  calibration  and 
discusses  the  development  and  operation  of  the  calibration 
program.  Chapter  V  presents  the  results  of  the  wavelength 
calibration  over  all  five  spectral  regions  studied  thus 
far.  Chapter  VI  presents  conclusions  and  recommends  areas 
for  future  research.  Finally,  the  relevant  software  code  is 
included  in  the  Appendix. 


2 


II.  BACKGROUND 


A.  SPECTRAL  IMAGING 

Spectral  imaging  is  a  rapidly  advancing  area  in  the 
field  of  remote  sensing  (Wolfe,  3)  .  Imaging  spectrometers 
have  combined  the  venerable  techniques  of  imaging  and 
spectroscopy,  thus  making  available  extremely  useful, 
integrated  data  sets  that  could  not  be  acquired  by  using 
either  technique  alone. 

Traditional  imaging  is  concerned  with  the  accurate 
measurement  of  light  intensity  over  a  two-dimensional 
scene.  Spatial  variations  in  intensity  are  used  to  detect 
features  and  patterns.  This  is  very  useful  for  the 

detection  of  objects,  but  there  are  certain  inherent 

limitations.  Spatial  features  such  as  size,  shape,  color, 
or  pattern  must  be  relied  upon  to  characterize  such 
objects.  Physical  attributes  such  as  material  composition 
can  be  difficult  to  determine  using  purely  image-based 
means.  Furthermore,  camouflage  can  be  used  with  varying 
degrees  of  success  to  interfere  with  detection  by  imaging 
sensors,  such  as  by  artificially  disguising  edge  features 
in  a  target  scene.  Decoys  with  properly  reproduced 

features  can  be  used  to  disrupt  the  accuracy  of  image 

interpretation.  These  drawbacks  can  therefore  result  in 
inaccurate  knowledge  of  an  imaged  scene. 

Spectroscopy  is  a  well-established  field  concerned 
with  the  study  of  variations  in  light  intensity  as  a 
function  of  wavelength  or  frequency.  Differing  materials 
exhibit  different  spectral  properties  due  to  their  atomic 
and  molecular  composition.  These  characteristics  include 
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material-specif ic  wavelengths  where  electromagnet ic  energy 
is  absorbed  (absorption  lines  or  bands)  or  emitted 
(emission  lines  or  bands) .  Under  favorable  circumstances, 
spectrometers  capable  of  detecting  such  spectral 
characteristics  can  be  used  to  determine  both  the  materials 
being  observed  and  some  characteristics  of  their 
environment,  such  as  ambient  temperature. 

Spectral  imaging  integrates  imaging  and  spectroscopy 
to  collect  voluminous  amounts  of  information  from  a  scene. 
Incident  light  from  a  scene  is  detected  and  recorded 
according  to  both  position  within  the  image  and  wavelength. 
A  spectral  image  is  thus  a  function  of  a  three-dimensional 
group  of  independent  variables,  typically  displayed  as  an 
abstract  hyperspectral  cube,  or  hypercube.  For  each 
spatial  element  (pixel)  of  the  image,  a  spectral  imager 
records  not  just  the  total  intensity  of  incident  light  in  a 
single  wavelength  band,  but  the  intensity  over  many  bands 
of  differing  wavelength.  The  hypercube  can  then  be 
dissected  and  analyzed  as  a  collection  of  two-dimensional 
images  at  various  wavelengths. 

Due  to  the  different  spectral  responses  of  various 
materials,  features  or  objects  previously  indistinguishable 
using  a  conventional  imager  can  be  detected  and  classified. 
Furthermore,  the  spectral  responses  of  individual  pixels 
can  be  examined  in  order  to  determine  the  material 
composition  and  environmental  conditions  in  those  small 
portions  of  an  image.  Careful  comparisons  with  nearby 
pixels  can  then  aid  in  the  identification  of  variations  in 
the  materials  present  in  an  emitting  target. 
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Figure  1.  Hyperspectral  Cube  of  Industrial  Smokestacks, 

Recorded  by  NUVIS  (From:  Marino, 10) 


Spectral  imagers  have  many  civilian  and  military 
applications  for  remote  sensing.  The  ability  to  determine 
different  material  compositions  across  a  scene  has 
tremendous  potential  for  Earth  sciences  and  resource 
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identification  purposes.  Likewise,  the  defeat  of  simple 
optical  camouflage  techniques  by  examining  many  regions  of 
the  electromagnetic  spectrum  is  a  valuable  capability  for 
military  sensors.  The  vast  potential  for  spectral  imaging 
devices  has  therefore  resulted  in  recent  research  into 
their  further  development.  One  product  of  this  research  is 
the  LINUS  project  at  the  Naval  Postgraduate  School. 

B .  LINUS 

The  LINUS  project  was  founded  following  the  successful 
use  of  its  predecessor,  NUVIS,  to  detect  and  measure  SO2 
concentrations  in  atmospheric  plumes.  The  NUVIS  optical 
design  was  limited  to  the  relatively  small  ultraviolet  (UV) 
spectral  band  of  approximately  300  to  375nm.  LINUS  was 
developed  to  apply  lessons  from  the  NUVIS  experience  and  to 
incorporate  several  enhancements,  such  as  the  instrument's 
tunability  across  the  near  UV  spectrum. 

In  order  to  detect  and  measure  many  other  potentially 
harmful  chemical  or  biological  agents,  LINUS  was  designed 
to  operate  with  greater  resolution  over  the  much  larger 
ultraviolet  spectral  range  from  200  to  400  nm.  For  the 
purposes  of  this  thesis  research,  this  region  has  been 
subdivided  into  five  smaller  bands  that  are  segregated  by 
the  use  of  individual  UV  filters. 

All  LINUS  components  have  been  selected  to  enable 
rugged  field  use  of  the  instrument.  The  actual  detector 
assembly  is  mounted  on  a  heavy-duty  tripod  and  sealed 
within  a  black  Acrylonitryl  Butadiene  Styrene  (ABS)  plastic 
shell  to  shield  the  optics  from  external  UV  sources  and 
environmental  contamination.  A  f ield-ruggedized  PC-type 
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data  acquisition  and  control  computer  and  power 
conditioning  hardware  have  been  mounted  in  shockproof 
transportation  containers  to  allow  their  use  in  the  field. 
While  all  developmental  work  to  date  has  been  conducted  in 
a  laboratory  environment,  LINUS  will  be  capable  of 
operation  in  harsher  outdoor  conditions. 


Figure  2 . 


LINUS  Field-Support  Equipment 


The  LINUS  optical  design  uses  an  image-intensified  CCD 
camera,  optimized  for  the  ultraviolet  region,  several 
compound  lenses,  a  planar  image-scanning  mirror,  and  a 
diffraction  grating  in  order  to  generate  spectral  image 
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data.  Incident  light  from  a  distant  target  is  reflected 
off  the  scanning  mirror  and  through  a  UV  bandpass  filter 
with  a  moderate  bandwidth  of  approximately  20  nanometers 
(nm)  .  The  UV  light  transmitted  through  the  filter  passes 
through  an  initial  objective  lens  system,  which  focuses  an 
image  of  the  distant  scene  onto  a  narrow  slit.  This  slit 
allows  only  a  thin  vertical  slice  of  the  scene  to  enter  the 
spectrometer.  During  data  recording,  successive  slices  are 
selected  by  rotating  the  scanning  mirror  under  computer 
control.  The  vertical  image  slice  then  passes  through  a 
collimating  lens  system  and  on  to  a  plane  diffraction 
grating.  The  grating  disperses  the  light  according  to  the 
grating  equation 

m^  =  d(sin  0i-sin0o)  [Eg.  2.1] 


where : 

m  =  Order  of  the  diffraction  set  (±1,  ±2,  ±3...), 
X  =  Wavelength  of  the  diffracted  light, 

d  =  Diffraction  grating  inter-ruling  spacing, 

0i  =  incident  angle, 

0o  =  output  angle . 


Since  m,  d,  and  0i  are  fixed,  light  of  different 
wavelengths  (X)  will  be  diffracted  through  different  angles 
(0o) .  The  result  is  a  horizontal  dispersion  of  the  incident 
light  corresponding  to  its  wavelength. 


The  light  is  then  focused  by  a  camera  objective  lens 
system  onto  a  UV-sensitive  microchannel  plate  coupled  to  a 
512x512  pixel  charge-coupling  device  (CCD)  detector  array. 


Because  the  optical  design  is  stigmatic,  the  vertical 
positions  of  light  incident  on  the  CCD  array  correspond  to 
the  vertical  location  of  corresponding  portions  of  the  thin 
slice  of  scene  being  imaged.  Thus  the  vertical  coordinates 
of  the  recorded  image's  pixels  are  in  one-to-one 
correspondence  with  the  vertical  coordinates  of  the  scene 
being  focused  on  the  slit.  The  horizontal  pixel  coordinate 
of  the  measured  image  corresponds  to  the  wavelength  of 
light  diffracted  from  each  vertical  element  of  the  scene. 
This  image  is  recorded  as  a  single  two-dimensional  slice  of 
the  final  hyperspectral  cube.  In  order  to  generate  the 
second,  horizontal  scene  dimension  required  for  a  full  2-D 
image,  the  mirror  is  slightly  rotated  by  a  precision 
closed-loop  DC  servo  system.  This  allows  an  adjacent 
vertical  slice  of  scene  to  be  imaged  by  the  camera.  The 
entire  spectral  image  is  built  up  by  interleaving  image 
frames  with  the  small  rotations  of  the  camera  mirror  until 
the  entire  scene  has  been  panned  across  the  camera' s 
instrument's  entrance  slit.  The  data  acquisition  computer 
then  stores  the  resultant  spectral  image  for  later 
analysis . 

If  the  instrument  is  operated  at  its  full  data 
resolution  of  512  horizontal  image  samples  by  512  vertical 
image  samples  by  512  wavelength  samples  by  12  bits  (2 
bytes)  per  pixel,  the  total  data  storage  requirement  for 
one  scene  will  be  512  x  512  x  212  x  2  =  268,435,456  bytes. 
In  practice  one  will  seldom  need  such  a  large  volume  of 
spectral  imaging  data. 
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LINUS 


Figure  3.  LINUS  Optical  Design 


10 


III.  EXPERIMENTAL  MEASUREMENTS 


A.  PURPOSE 

In  order  for  the  information  from  a  hypercube  of 
spectral  data  to  be  of  maximum  usefulness,  the  dimensions 
of  all  three  axes  must  be  accurately  calibrated.  The  two 
physical  image  dimensions  (x  and  y)  can  be  calibrated  from 
the  properties  (i.e.  effective  focal  lengths)  of  the 
optical  components  and  their  relative  spacing  within  the 
instrument.  This  spatial  calibration  can  be  verified  by 
imaging  an  object  of  known  size  and  distance  and 

determining  its  resultant  size  in  terms  of  pixel  width  and 
height.  The  third,  wavelength  (X),  axis  can  also  be 
calibrated  in  principle  by  using  the  grating  equation  and 
the  focal  properties  of  the  optical  train.  In  practice, 
however,  such  calculations  are  useful  only  for  rough 
approximations  of  measured  wavelengths.  There  are  two 
primary  reasons  for  this.  First,  the  use  of  chromatic 
lenses  in  the  optics  causes  variations  in  focal  length,  and 
hence  in  image  scale,  because  of  chromatic  aberration 
(Pedrotti  and  Pedrotti,  102)  .  The  large  spectral  range 
over  which  LINUS  is  designed  to  operate  aggravates  this 
effect.  Second,  the  inherent  non-linearity  of  the 

diffraction  grating  equation  (Eqn.  2.1)  leads  to  a 

nonlinear  wavelength  scale  across  the  recorded  image  plane. 
In  order  to  perform  spectral  absorption  analysis  and  other 
meaningful  wavelength-based  studies,  wavelength  must  be 
precisely  known.  This  thesis  project  was  therefore 

initiated  to  develop  a  method  for  wavelength  calibration  by 
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using  the  LINUS  optical  train  as  a  spectrometer  to  observe 
an  independently  calibrated  spectral  standard  source. 

B .  EXPERIMENTAL  SETUP 

In  order  to  conduct  this  wavelength  calibration 
research,  LINUS' s  ability  to  perform  two-dimensional 
imagery  was  not  needed.  The  instrument  was  configured  in  a 
purely  spectroscopic  mode.  Readers  who  are  interested  in 
LINUS'S  image  scanning  subsystem  can  find  a  thorough 
description  in  the  earlier  thesis  by  R.  Kompatzki 
(Kompatzki,  25)  .  For  convenient  access  during  the 
experiments,  all  of  LINUS' s  optical  and  optomechanical 
components,  except  the  image  scanner,  were  removed  and 
positioned  on  an  optical  table. 

A  platinum  (Pt)  hollow  cathode  lamp  was  used  as  the 
spectral  calibration  standard  source.  Atomic  platinum 
exhibits  a  rich  spectrum  of  emission  lines  across  the 
entire  near-ultraviolet  region.  The  National  Institute  of 
Standards  and  Technology  (NIST)  has  catalogued  the  precise 
wavelengths  and  relative  intensities  of  these  emission 
lines,  and  has  made  the  data  readily  available  on  their 
website  (Sansonetti) .  This  has  made  platinum  hollow 
cathode  lamps  one  of  the  most  commonly  used  wavelength 
standards  for  ultraviolet  analysis.  Additionally,  a 
regulated  DC  power  supply  was  built  to  consistently  deliver 
the  optimum  rated  current  to  the  platinum  lamp.  As  shown 
in  Figure  5,  the  lamp  was  placed  along  the  optical  axis  of 
the  instrument,  in  approximately  the  position  that  would 
normally  be  occupied  by  the  scanning  mirror. 
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Figure  4.  Pt  Hollow-Cathode  Lamp  and  DC  Power  Supply 


Figure  5.  Placement  of  Pt  Lamp  Along  Optical  Axis 
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Once  the  LINUS  optical  components  and  platinum  lamp 


were  mounted  to  the  table,  simple  light-blocking  baffles 
were  installed  to  reduce  the  amount  of  stray  ultraviolet 
light  reaching  the  camera  to  an  acceptable  level.  The 
different  UV  spectral  regions  were  selected  by  using 
filters  with  different  pass  bands.  The  filter 

characteristics  are  summarized  in  Table  1,  and  the 
manufacturer' s  transmission  curves  are  shown  in  Figures  6 
through  10  (Omega  Optical  Company) . 


Filter  Number 

Manufacturer 

Peak 

Passband 

Designation 

Transmission 

(50%  of  Peak) 

Wavelength (nm) 

(nm) 

1 

220BP10 

220 

216-226 

2 

289BP10 

286 

283-293 

3 

300BP10 

296 

293-304 

4 

334BP10 

330 

326-346 

5 

370BP10 

367 

363-375 

Table  1.  UV  Filter  Characteristics 
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Figure  7  . 


Filter  Transmission  (%T)  Curve:  289BP10 


Figure 
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Figure  9. 


Filter  Transmission  (%T)  Curve:  334BP10 
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Figure  10.  Filter  Transmission  (%T)  Curve:  370BP10 

C .  DATA  COLLECTION 

Before  any  data  was  taken,  several  exercises  were 
conducted  for  familiarization  with  the  data  acquisition 
software.  First,  one  of  the  UV  bandpass  filters  was 
selected  and  installed.  Then,  after  the  proper  operation 
of  the  equipment  was  verified,  the  platinum  lamp  was 
energized  and  raised  to  its  specified  operating  current 
level.  The  diffraction  grating  was  rotated  so  that  the 
zeroth  order  image  of  the  slit  was  properly  focused  onto 
the  camera's  image  intensifier.  The  slit  width  was 
narrowed  until  its  image  appeared  crisp  and  the  measured 
intensity  level  at  the  camera  was  acceptable,  without 
either  excessive  noise  or  intensity  saturation. 

The  diffraction  grating  was  then  rotated  until  the 
first  order  set  of  diffraction  lines  became  centered  in  the 
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camera's  field  of  view.  Rough  focusing  was  accomplished  by 
alternately  adjusting  the  positions  of  the  collimating  and 
objective  lenses  with  respect  to  the  slit  and  the  camera 
focal  plane,  respectively.  Finer  adjustment  of  the  focus 
was  accomplished  by  using  the  camera' s  image  acquisition 
software,  WinView/32  (Roper  Scientific),  to  enlarge  a  small 
region  of  the  viewable  area  containing  several  weak 
emission  lines  and  further  refining  the  lens'  positions  to 
sharpen  the  image . 

Once  the  diffracted  emission  lines  were  properly 
focused,  the  image  was  recorded  and  saved  in  Tagged  Image 
File  (.TIF)  format.  This  image  was  subsequently  exported 
to  and  processed  by  the  calibration  program  written  with 
the  MATLAB  software  package.  The  procedure  was  then 
repeated  for  the  remaining  UV  filters.  Representative 
images  for  each  of  the  five  filters  are  shown  in  Figures  11 
through  15. 

In  these  figures,  as  discussed  earlier,  the  vertical 
dimension  corresponds  to  vertical  spatial  information, 
while  the  horizontal  dimension  is  the  spectral  dispersion 
direction . 
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Figure  1 3 . 


Figure  14 . 
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IV.  DATA  ANALYSIS 


A.  OVERVIEW 

The  bulk  of  the  work  involved  in  this  thesis  project 
consisted  of  creating  a  MATLAB  software  program  capable  of 
performing  pattern-matching  operations  between  measured 
spectra  and  corresponding  known  or  standard  spectra.  The 
desired  product  was  a  calibration  curve  that  relates  image 
pixel  number  versus  wavelength  across  each  observed 
spectral  band.  Cross-correlation,  which  is  a  familiar 
mathematical  method  for  comparing  the  similarity  of  two 
different  patterns,  was  used  to  optimize  the  alignment 
between  the  measured  and  standard  spectral  data.  Spectral 
line  data  recorded  by  LINUS  were  compared  to  expected 
values  from  the  tabulated  NIST  Pt  standard  line  list 
(Sansonetti) .  The  NIST  wavelength  scale  was  incrementally 
contracted  or  elongated,  then  shifted  until  the  LINUS  and 
standard  emission  lines  were  optimally  matched,  essentially 
performing  an  affine  transformation  on  the  wavelength 
scale.  Once  the  Pt  standard  spectral  lines,  for  which 
relative  intensity  and  absolute  wavelength  were  known,  were 
properly  aligned  with  the  LINUS  spectral  lines  the  known 
wavelengths  could  be  associated  with  the  LINUS  lines. 
Finally,  relating  the  pixel  location  of  the  spectral  lines 
to  their  corresponding  wavelengths  completed  the  wavelength 
calibration  process. 

B .  COMPUTATIONAL  THEORY 

1.  Preprocessing  of  Platinum  Standard  Spectral  Data 

Significant  preprocessing  of  the  NIST  spectral  data 
was  required  in  order  for  a  valid  correlation  to  be  made 
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with  the  measured  LINUS  data.  This  preprocessing  included 
elimination  of  irrelevant  data  and  simulating  the 
attenuation  of  the  UV  bandpass  filters  within  each  spectral 
band  that  was  studied. 

The  tabulated  data  contained  wavelength  and  relative 
intensity  information  on  over  3200  emission  lines  between 
the  wavelengths  of  2000  and  4000  angstroms.  Given  the 
large  calculat ional  requirements  of  repeatedly  conducting 
cross-correlation  within  iteration  loops,  inclusion  of 
irrelevant  spectral  lines  would  have  been  extremely 
inefficient  and  computationally  slow.  Therefore,  after 
loading  the  NIST  data  from  a  ".dat"  data  file,  all  values 
significantly  outside  of  the  given  UV  filter  passband  were 
discarded.  Furthermore,  the  user  could  specify  a  minimum 
intensity  threshold  below  which  weaker  lines  that  were  less 
discernable  by  LINUS  could  be  discarded.  The  resultant 
data  filtering  typically  reduced  the  number  of  lines 
considered  during  a  given  calibration  to  between  50  and 
100,  which  was  much  less  processor-intensive. 

Since  the  tabulated  NIST  data  do  not  include  actual 
line  profiles,  but  only  peak  intensities,  a  Gaussian  line 
shape  function  was  used  to  broaden  each  of  the  tabulated 
standard  spectral  lines  via  simple  numeric  convolution. 
Neighboring  lines,  particularly  those  whose  wavelengths 
differed  by  an  amount  less  than  or  equal  to  the 
characteristic  width  of  the  simulated  line  shape  function, 
were  blended  together  by  this  process.  This  blending  was 
also  expected  in  the  experimental  LINUS  spectrum. 

Additional  preprocessing  was  necessary  to  account  for 
the  transmissibility  effects  of  the  UV  filters. 
Transmissibility  versus  wavelength  was  included  for  many 


data  points  within  each  filters'  passband.  Interpolation 
for  each  relevant  NIST  emission  line  determined  the 
expected  transmissibility  at  that  wavelength.  The  relative 
intensity  of  each  line  was  then  modified  by  its 
corresponding  transmission  factor.  The  resultant  data 
represented  the  expected  relative  intensities  of  spectral 
lines  measured  by  LINUS. 
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Figure  16.  Removal  of  Low  Intensity  Emission  Lines 
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Figure  17.  Simulated  Broadening  and  Filtering  of 

Representative  Platinum  Spectra  in  the  3200-3700  Angstrom 

Spectral  Region 


2 .  Preprocessing  of  LINUS  Data 


The  filtered 

and 

spectrally 

dispersed 

light  from 

the 

laboratory 

platinum  hollow 

cathode  lamp 

formed  a 

two- 

dimensional 

image 

on 

the  LINUS 

camera  focal  plane. 

In 

order  to  compare 

data 

from 

this 

image  with 

intensity 

and 

wavelength 

data 

from 

the 

NIST 

standard. 

the  vertical 

spatial  information  in  each  image  was  not  needed. 
Therefore,  to  improve  signal-to-noise  characteristics  of 
the  spectrum,  the  light  intensity  reaching  each  element  in 
a  vertical  pixel  row  was  summed,  or  coadded.  In  addition, 
the  user  was  able  to  specify  an  intensity  threshold  in 
order  to  account  for  the  intensity-independent  noise  floor 
across  the  CCD  array.  The  result  was  a  one-dimensional 
measured  spectral  intensity  value  versus  image  pixel  row 
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address,  stored  in  a  data  array  that  was  suitable  for 
subsequent  correlation  with  the  corresponding  array  of  NIST 
standard  spectral  data  as  outlined  above. 


Figure  18.  LINUS  Image  Capture  and  Preprocessing: 

Top  Frame:  Raw  2-Dimensional  Data 
Middle  Frame:  Noise  Threshold  Definition 
Bottom  Frame:  Resultant  Preprocessed  Data 

3.  Pattern  Matching  and  Cross-correlation 

Cross-correlation  is  a  method  of  measuring  the  degree 
of  similarity  between  different  signals.  In  the  context  of 
this  thesis  project,  two  conditions  must  be  met  in  order 
for  two  signals  to  reach  maximum  agreement.  First,  the 
horizontal  data  scale,  such  as  time  or  frequency,  must  be 
equal.  For  this  project  the  horizontal  scale  was 

wavelength,  scaled  as  described  above.  Second,  the  two 
signals  must  be  optimally  aligned  with  each  other,  i.e.  any 
lateral  shifts  between  the  emission  lines  must  be 
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minimized.  When  both  of  these  conditions  are  met,  the 
cross-correlation  of  two  signals  will  reach  a  maximum 
value.  The  equation  for  cross-correlation  as  a  function  of 
data  series  offset,  r (d) ,  is  given  by:  (Bourke) 

r(d)  =  -j^‘=^  — I  ^  [Eg.  3.1] 

*J^[(y(i-d)-yf 

where 

x(i)  =  first  data  series  being  compared, 

y(i)  =  second  data  series  being  compared, 

i  =  index  of  the  data  sample  (1,2,3...  n)  , 

n  =  length  of  the  data  samples, 

d  =  sample  offset  between  x(i)  and  y(i) 

(0,l,...n-l)  . 

The  second  signal  is  shifted  relative  to  the  first  by 
iterating  the  delay  term,  d,  yielding  a  maximum  value  of 
cross-correlation  when  the  signals  are  best  aligned. 

The  MATLAB  one-dimensional  cross-correlation  function 
"xcorrl"  was  implemented  to  compare  the  preprocessed  LINUS 
data  with  the  modified  NIST  spectrum.  Because  the 

wavelength  scale  of  the  LINUS  data  was  unknown  a  priori,  an 
algorithm  was  developed  to  perform  an  affine  transform  on 
the  NIST  data  wavelength  scale  by  iteration  of  the 
horizontal  scale  through  a  specified  band.  At  each 
iteration  of  wavelength  scaling,  the  peak  cross-correlation 
value  was  recorded.  Plotting  the  maximum  cross-correlation 
value  against  the  scale  elongation  factor  enabled  the  point 
for  which  the  LINUS  and  NIST  data  best  agreed  to  be  readily 
identified. 
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Figure  19.  Cross-Correlation  Peak  vs.  Elongation  Factor 
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Figure  20.  Correlogram  Showing  Peak  at  Best  Shift 
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This  "best  scale  factor"  was  used  to  reiterate  the 
cross-correlation  and  to  determine  the  optimum  shift 
required  to  maximize  the  cross-correlation.  Application  of 
this  shift  to  the  horizontal  axis  of  LINUS  data,  i.e.  pixel 
number,  finally  resulted  in  the  alignment  of  NIST  and  LINUS 
spectral  lines  with  equally  scaled  spectral  wavelengths. 
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Figure  21.  Aligned  LINUS  and  NIST  Data 


The  LINUS  horizontal  scale, 
therefore  be  related  to  the 
wavelength . 


pixel  row  number,  could 
NIST  horizontal  scale. 
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Figure  22.  Wavelength  Calibration  Curve 


The  reader  must  be  cautioned  that  such  a  linear 
calibration  will  necessarily  be  valid  over  only  a  limited 
range  of  wavelengths,  and  that  it  must  be  performed  for 
each  spectral  region.  The  diffraction  grating  equation  is 
nonlinear  (sinusoidal) ,  meaning  that  the  wavelength  scale 
will  also  be  nonlinear. 

A  widely  accepted  measure  of  goodness-of-f it  between 

2 

two  distributions  is  the  %  parameter.  (Taylor)  Once  the 

cross-correlated  LINUS  and  standard  spectra  have  been 
2 

matched,  a  %  value  may  be  calculated  to  characterize  the 

2 

agreement.  The  %  parameter  is  generated  as  follows: 
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[Eqn .  4.1] 


1=1 

where  ai  and  bi  are  the  corresponding  samples  of  the 
distributions,  and  the  summation  is  over  all  relevant 
samples.  For  the  purpose  of  this  research,  ai  was  the  LINUS 
spectrum  of  interest,  and  bi  was  the  corresponding  NIST 
spectrum . 

2 

Analysis  of  the  %  parameter  was  conducted  to  verify 

the  pattern-matching  capability  of  the  cross-correlation 

program.  A  "best  fit"  between  two  distributions  exists 
2 

when  the  %  parameter  is  minimized.  This  value  was  first 
calculated  for  the  matched  intensity  distributions.  The 
alignment  of  the  intensity  distributions  was  then  shifted 

2 

by  user-specified  amounts.  The  %  parameter  was  calculated 

for  the  manually  offset  intensity  distributions  and 

2 

verified  to  be  greater  than  the  %  value  for  the  aligned 
distributions . 
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V.  RESULTS 


A.  OVERVIEW 

Separate  wavelength  calibrations  were  performed  for 
the  ultraviolet  bands  associated  with  each  of  the  five 
filters.  In  each  case,  the  automated  results  were 

confirmed  for  proper  alignment  with  the  NIST  spectral  data 
by  visual  comparison  and  by  analysis  of  the  parameter. 
The  minimum  parameters  for  each  spectral  band  are 

presented  in  tabular  form. 


Filter 

Minimum  Value 

220nm  Filter 

6.16  E7 

289nm  Filter 

9.59  E6 

300nm  Filter 

1.29  E5 

334nm  Filter 

1.63  E5 

370nm  Filter 

1.44  E3 

Table  2.  Minimum  Parameters  for  Calibration  Results 

Associated  with  Spectral  Region  of  Each  Filter 


B.  220NM  FILTER 

In  this  wavelength  range,  the  512  pixel  rows  of  the 
camera' s  focal  plane  were  found  to  span  a  wavelength 
interval  of  51.1  nm,  from  197.8  to  248. 9nm. 
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Figure  23.  Pt  Standard  Preprocessing:  220nm  Filter 

Top  Panel:  Pt  Standard  Emission  Lines 
Middle  Panel:  220  nm  Filter  Curve 
Bottom  Panel:  Final  Preprocessed  Output 


Raw  LINUS  Image  of  Pt  Spectrum  Using  220  nm 
Filter  (Axes  Correspond  to  Pixel  Numbers) 
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Figure  24 . 
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Figure  25.  Coadded  Pt  Spectrum  and  Threshold  Subtraction 

Using  220nm  Filter 
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Figure  26.  Data  Alignment  Using  220nm  Filter 

Top  Panel:  Preprocessed  NIST  Reference  Data 
Bottom  Panel:  Aligned  LINUS  Data 
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Figure  2 7 . 


Wavelength  Calibration  Curve  Using  220  nm  Filter 


C.  289NM  FILTER 


The  camera' s 
wavelength  interval 


focal 
of  22.0 


plane  was  found  to 
nm,  from  274.1  to  296.1 


span 
nm . 
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Figure  28.  Pt  Standard  Preprocessing:  289nm  Filter 

Top  Panel:  Pt  Standard  Emission  Lines 
Middle  Panel:  289  nm  Filter  Curve 
Bottom  Panel:  Final  Preprocessed  Output 


Raw  LINUS  Image  of  Pt  Spectrum  Using  289  nm 
Filter  (Axes  Correspond  to  Pixel  Numbers) 
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Figure  29 


Figure  30.  Coadded  Pt  Spectrum  and  Threshold  Subtraction 

Using  289nm  Filter 
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Figure  31.  Data  Alignment  Using  289nm  Filter 

Top  Panel:  Preprocessed  NIST  Reference  Data 
Bottom  Panel:  Aligned  LINUS  Data 
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Figure  32.  Wavelength  Calibration  Curve  Using  289  nm  Filter 
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D.  300NM  FILTER 

The  camera' s  focal  plane  was  found  to  span  a 
wavelength  interval  of  34.5  nm,  from  282.4  to  316.9  nm. 
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Figure  33.  Pt  Standard  Preprocessing:  300nm  Filter 

Top  Panel:  Pt  Standard  Emission  Lines 
Middle  Panel:  300  nm  Filter  Curve 
Bottom  Panel:  Final  Preprocessed  Output 


Raw  LINUS  Image  of  Pt  Spectrum  Using  300  nm 
Filter  (Axes  Correspond  to  Pixel  Numbers) 
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Figure  34 . 
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Figure  35.  Coadded  Pt  Spectrum  and  Threshold  Subtraction 

Using  300nm  Filter 
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Figure  36.  Data  Alignment  Using  300nm  Filter 

Top  Panel:  Preprocessed  NIST  Reference  Data 
Bottom  Panel:  Aligned  LINUS  Data 
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Figure  37.  Wavelength  Calibration  Curve  Using  300  nm  Filter 


E.  334NM  FILTER 


The  camera' s 
wavelength  interval 


focal 
of  34.5 


plane  was  found  to 
nm,  from  326.2  to  360.7 


span 

nm. 
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Figure 


Figure  39 


38.  Pt  Standard  Preprocessing:  334nm  Filter 

Top  Panel:  Pt  Standard  Emission  Lines 
Middle  Panel:  334  nm  Filter  Curve 
Bottom  Panel:  Final  Preprocessed  Output 


Raw  LINUS  Image  of  Pt  Spectrum  Using  334  nm 
Filter  (Axes  Correspond  to  Pixel  Numbers) 
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Figure  40.  Coadded  Pt  Spectrum  and  Threshold  Subtraction 

Using  334nm  Filter 


Figure  41.  Data  Alignment  Using  334nm  Filter 

Top  Panel:  Preprocessed  NIST  Reference  Data 
Bottom  Panel:  Aligned  LINUS  Data 
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Figure  42.  Wavelength  Calibration  Curve  using  334  nm  Filter 


F.  370NM  FILTER 


The  camera' s 
wavelength  interval 


focal 
of  34.1 


plane  was  found  to  span 
nm,  from  345.8  to  379.9  nm. 
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Figure 


Figure  44 


43.  Pt  Standard  Preprocessing:  370nm  Filter 

Top  Panel:  Pt  Standard  Emission  Lines 
Middle  Panel:  370  nm  Filter  Curve 
Bottom  Panel:  Final  Preprocessed  Output 


Raw  LINUS  Image  of  Pt  Spectrum  Using  370  nm 
Filter  (Axes  Correspond  to  Pixel  Numbers) 

46 


f  10^  irt<*i9njr»T 


Figure  45.  Coadded  Pt  Spectrum  and  Threshold  Subtraction 

Using  370nm  Filter 
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Figure  46.  Data  Alignment  Using  370nm  Filter 

Top  Panel:  Preprocessed  NIST  Reference  Data 
Bottom  Panel:  Aligned  LINUS  Data 
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Figure  47 . 


Wavelength  Calibration  Curve  Using  370  nm  Filter 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  project  was  successful  in  generating  wavelength 
calibration  data  for  the  LINUS  instrument.  The  automatic 
wavelength  calibration  for  all  five  spectral  bands  was 
easily  verified  qualitatively  by  eye.  Quantitatively,  the 
aligned  observed  LINUS  data  generated  minimal  error 

values  when  compared  to  the  platinum  standard. 

The  calibration  program  as  completed  is  capable  of 
generating  accurate  calibration  curves  for  all  of  the 
spectral  bands  studied  to  date.  These  five  spectral  bands 
do  not,  however,  completely  span  the  entire  LINUS  design 
operating  band  of  200-400nm.  It  is  therefore  recommended 
that  additional  filters  be  procured  to  span  these  gaps  in 
the  LINUS  operating  band.  Once  the  calibration  procedure 
has  been  conducted  for  the  remaining  filters,  LINUS  will  be 
capable  of  detecting  and  quantifying  a  large  number  of 
chemical  species. 
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APPENDIX:  WAVELENGTH  CALIBRATION  PROGRAM  CODE 


0000000000000000000000000000000000000000000000000000 


LINUS  Wavelength  Calibration  Program 


2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2-2- 

0000000000000000000000000000000000000000000000000000 


This  program  integrates  previously  developed 
subprograms  in  order  to  perform  wavelength 
calibration  of  the  LINUS  UV  imaging  spectrometer. 

The  process  consists  of: 

1.  Preprocessing  NIST  data 

Loading  from  . dat  file 
Pre-filtering  for  intensity,  filter 


passband 

%  Correcting  for  UV  filter  transmission 

%  Broadening  the  spectral  lines 

%  Reindexing  for  x-axis  graduation 

o, 

o 

%  2.  Preprocessing  LINUS  data 

%  Loading  from  .tif  file 

%  Converting  2-D  image  to  1-D  intensity 

%  Pre-filtering  for  intensity 

o, 

o 

%  3.  Cross-correlation 

%  "Stretch"  NIST  for  maximum  xcorr  value 

%  Using  "Best  stretch",  find  best  shift 

between  NIST  and  LINUS 

%  Complete  calibration  -  wavelength  vs.  pixel 

line  number 


VERSION  NOTES: 


clear 

%  USER  SPECIEIED  PARAMETERS 

filterChoice  =  1%  l=220nm,  2=289nm,  3=300nm,  4=334nm, 

5=37  Onm 

flipLINUS  =0  %  0  if  image  normal,  1  if  image 

flipped 

%  289nm  data  only  (to  date) 

if  ( f ilterChoice==l )  %  220nm 

MinWavelength  =  2100; 

% [Angstroms ] 
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2400; 


MaxWave length 
% [Angstroms ] 

IntensityThreshold  =  2.5e4;  %3e4 

LINUSIntensityThreshold  =  7e4  %2e4 

LINUSfile  =  'a220nmlst.tif' 

%  Input  filter  transmission  data 

y  =  [0.00  0.010  0.04  0.13  0.320  0.33  0.300  0.210  0.040 

0.010  0.000]'; 

t  =  10*[210  212  215  217  220  221  223  225  230  235  240]'; 


%  Affine  Transformation  Range 
Init lalStretchFactor  =  . 14% . 20% . 25; % . 14 
FinalStretchFactor  =  . 17% . 25% . 2 9; % . 1 6 
Stretchinterval  =  .005; 

end 


if  (f ilterChoice==2)  %  289nm 

MinWavelength  =  2700; 

% [Angstroms ] 

MaxWavelength  =  3100; 

% [Angstroms ] 

IntensityThreshold  =  2e5; 

LINUSIntensityThreshold  =  4e4; 

LINUSfile  =  'a289nmlstMoreMCP.TIF' 
flipLINUS  =1  %  0  if  image  normal 

flipped 


%8e4 

%4e4 

1  if  image 


%  289nm  data  only  (to  date) 
Input  filter  transmission  data 


y  =  [0.000  0.010  0.020  0.070  0.300  0.370  0.270  0.120 

0.010  0.000  0.000]'; 

t  =  10*[270  272  275  280  285  288  290  295  300  305  310]'; 


%  Affine  Transformation  Range 
Init lalStretchFactor  =  . 4 63% . 4 6; % . 50 
FinalStretchFactor  =  .5%. 50;%. 52 
Stretchinterval  =  .001; 

end 


if  (f ilterChoice==3)  %  300nm 

MinWavelength  =  2850; 

% [Angstroms ] 

MaxWavelength  =  3100; 

% [Angstroms ] 

IntensityThreshold  =  2e4; 

LINUSIntensityThreshold  =  5e4; 


LINUSfile  =  'a300nmlstFoc2.tif' 

%  Input  filter  transmission  data 


%3e4 

%6e4 
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y  =  [0.000  0.001  0.050  0.220  0.320  0.280  0.110  0.000]'; 
t  =  10*[280  285  290  295  297  300  305  310]'; 

%  Affine  Transformation  Range 
Init ialStretchFactor  =  .12;%. 12 
FinalStretchFactor  =  .19;%. 14 
Stretchinterval  =  .005; 

end 

if  ( f ilterChoice==4 )  %  334nm 

MinWavelength  =  3200; 

% [Angstroms ] 

MaxWavelength  =  3700; 

% [Angstroms ] 

IntensityThreshold  =  5e4;  %3e4 

LINUSIntensityThreshold  =  3e4;  %2e4 

LINUSfile  =  'a334nmlst.tif' 

%  Input  filter  transmission  data 
y  =  [0.000  0.003  0.140  0.450  0.330 

0.300  0.250  0.130  0.060  0.020 

0.000  0 ] ' ; 

t  =  10*[315  320  325  330  335  340  345  350  355  360  365  370]'; 

%  Affine  Transformation  Range 
Init IalStretchFactor  =  .3;%.37BSF 
FinalStretchFactor  =  .7;% 

Stretchinterval  =  .005; 

end 

if  ( f ilterChoice==5 )  %370nm 

MinWavelength  =  3400; 

% [Angstroms ] 

MaxWavelength  =  3900; 

% [Angstroms ] 

IntensityThreshold  =  5e4;  %3e4 

LINUSIntensityThreshold  =  6e4;  %2e4 

LINUSfile  =  'a370nmlst.tif' 

%  Input  filter  transmission  data 

y  =  [0  .003  .01  .03  .07  .13  .2  .25  .30  .32  .30  .28  .26 

.23  .20  .16  .12  .05  .03  .01  0  0]  '  ; 

t  =  1000*[3.40  3.45  3.50  3.55  3.60  3.62  3.63  3.64  3.65 
3.67  3.69  3.70  3.71  3.72  3.73  3.74  3.75  3.80  3.85  3.90  3.95 
4.00]  '  ; 

%  Affine  Transformation  Range 
Init IalStretchFactor  =  .35;%. 37 
FinalStretchFactor  =  .41;%. 29 
Stretchinterval  =  .005; 
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end 


%  Other  Definitions: 

%  Sampling  parameter  for  NIST  elongation 
divisions  =  2000; 

%  define  gaussian  curve  shape 

gauss  =  [.002  .018  .105  .368  .779  1  .779  .368 

.018  .002]; 


00000000000000000000000000000000000000000000000000000 

g,  g,  o,  g,  o,  g, 

0  0  0  0  0  0 


%  1.1  Data  import  and  assignment  into  two  Nxl  matri 

load  NISTFull.dat; 

NISTWavelength  =  NISTFull  (  :  ,  1 ) ; 

NISTIntensity  =  NISTFull  (:, 2 ) ; 

NISTlength  =  length (NISTWavelength) 
waveband  =  MinWavelength : 1 : MaxWavelength; 

figure  ( 1 ) 
subplot  (4,1,1)  ; 

plot (NISTWavelength, NISTIntensity) 
hold  on 

xdatum=MinWavelength : 1 : MaxWavelength; 
plot (xdatum, IntensityThreshold,  ' -r '  )  ; 
hold  off 

axis ( [MinWavelength, MaxWavelength, 0, 4e5] )  ; 


00000000000000000000000000000000000 
%%%  FINAL  DISPLAY  PLOT 
figure (40) 
subplot  (3,1,1); 

plot (NISTWavelength, NISTIntensity) 
hold  on 

xdatum=MinWavelength : 1 : MaxWavelength; 
plot (xdatum, IntensityThreshold,  ' -r ' ) ; 
hold  off 

axis ( [MinWavelength, MaxWavelength, 0, 4e5] ) %2e6 
xlabel ( 'Wavelength  (Angstroms) ' ) ; 
ylabel ( ' Intensity ' )  ; 

0000000000000000000000000000000000000000000 
hold  on 
figure  ( 1 ) 


.105 


g,  g,  g,  g,  g,  g, 
0  0  0  0  0  0 


ces 
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%  1.2  Discard  values  outside  of  bandwidth  and  below 

threshold 
1=1; 

while  i  <=  NISTlength 

if  (NISTWavelength ( i ) <MinWavelength  | 

NISTWavelength (i) >MaxWavelength  | 

NISTIntensity (i) <IntensityThreshold) 

NISTIntensity (!)=[] ; 

NISTWavelength (!)  =  []  ; 

NISTlength=NISTlength-l; 

i=i-l; 

end 
i=i+l ; 

end 

NISTlength  =  length (NISTWavelength) 
subplot  (4,1,2); 

plot (NISTWavelength, NISTIntensity) 

%  1.3  Gaussian  broaden  each  spectral  line 

NISTwl  =  NISTWavelength;  %  Irregularly  spaced, 
fractional  Angstroms 
NISTint  =  NISTIntensity; 

NISTwl  =  round (NISTwl ) ;  %  Rounds  to  integer  number 

of  Angstroms 


%f  igure  (2  0) 
subplot  (4,1,3) 
plot (gauss) 

axis ( [ 1 , length (gauss ) ,0,1]) 

newNISTwl  =  zeros (1, length (NISTwl) *length (gauss) ) ; 
newNISTint  =  zeros (1, length (NISTwl) *length (gauss) ) ; 
halfwidth  =  (length (gauss) -1) /2; 

A  =  1; 

while (A<=length (NISTint) ) 

B  =  1; 

while (B<=length (gauss) ) 

C= (A-1) * length (gauss) +B; 

newNISTwl (C)  =  NISTwl (A) -halfwidth+B; 

newNISTint (C) =  gauss (B) *NISTint (A) ; 

B=B+1; 

end 
A=A+ 1 ; 
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end 


newNISTwl;  %  Broadened,  overlapping  values 
newNISTint ; 

subplot  (4,1,4) 

plot (newNISTwl , newNISTint , ' r ' ) 

%  1.4  Sum  overlapping  intensity  distributions 

F  =[]; 

G=l; 

H  =  1; 

J  =  [  ]  ; 


while  (G<=length (waveband) ) %was  spectrum 
while  (H<=length (newNISTwl ) ) 

F  (H)  =  isequal (newNISTwl (H) , waveband (G) ) ; 
H=H+1; 

end 

J (G) =sum (newNISTint  (F) ) ; 

H=l; 

G=G+1; 

end 

J; 

figure  (2 ) 
subplot  (3,1,1) 

plot (newNISTwl , newNISTint , ' r ' ) 
hold  on 

plot  (waveband, J) 
hold  off 

NISTWavelength  =  waveband; 

NISTIntensity  =  J; 

%  1.3  NIST  SPECTRUM  MODIFICATION 

%  I.  3.  I  SPLINE  CURVE  FIT 

%  Interpolate  transmissability  at  each  relevant 

wavelength 

T2=NISTWavelength; 

FilterT  =  spline (t,y,T2); 


subplot  (3,1,2) ; 
plot (t, y, 'o') 

56 


axis ( [MinWavelength, MaxWavelength, 0 , . 5 ] ) ; 
hold  on 

plot(T2,FilterT, 'r. ') 

%  1.3.2  COMPUTE  EXPECTED  POST-EILTER  INTENSITY  VALUES 

Eilteredintensity  =  EilterT.*J; 
subplot  (3,1,3); 

plot (NISTWavelength, Eilteredintensity ,  'g'  ) 
axis ( [MinWavelength, MaxWavelength, 0, le5] )  ; 


00000000000000000000000000000000000 
%%%  FINAL  DISPLAY  PLOT 
figure (40) 
subplot  (3,1,2); 

%plot (t , y , 'o') 

axis ( [MinWavelength, MaxWavelength, 0 ,  . 5 ] )  ; 
hold  on 

plot  (T2,FilterT) 

xlabel  (  'Wavelength  (Angstroms)  ' )  ; 
ylabel  (  ' Transmissibility '  )  ; 

subplot  (3,1,3); 

plot  (NISTWavelength, Filteredintensity) 

axis ( [MinWavelength, MaxWavelength, 0, le5] ) ;%2e5, 3e5, le6, 4e5 
xlabel ( 'Wavelength  (Angstroms)  '  )  ; 
ylabel ( ' Intensity '  )  ; 

0000000000000000000000000000000000000000000 


%  2  PROCESS  LINUS  SPECTRAL  RESPONSE 

%  2.1  Data  import  and  assignment  into  2  Nxl  matrices 

LINUSimage  =  imread (LINUSf ile)  ; 
if (flipLINUS==l) 

LINUSimage  =  fliplr (LINUSimage) ; 

end 

figure  ( 3 ) 
subplot  (3,1,1); 
contour (LINUSimage) 


000000000000000 

o,  o,  o,  o,  o,  o, 

0  0  0  0  0  0 
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figure  (41) 
contour (LINUSimage) 
figure  ( 3 ) 

000000000000000000000 


LINUSintensity  =  sum (LINUSimage, 1 ) ; 
DATAlength  =  length (LINUSintensity) ; 

subplot  (3,1,2); 
hold  on 

plot  (LINUSintensity)  ; 
xxdatum=l : 1 : 512  ; 

plot (xxdatum, LINUSIntensityThreshold,  ' -r ' ) ; 
hold  off 

%Title (' Intensity  vs.  Pixel  row'); 
xlabel  (' Pixel  line  number'); 
ylabel ( ' Intensity '  )  ; 
axis  (  [ 0 , DATAlength, 0 , 5e5 ] )  ; 


000000000000000 

g,  o,  g,  o,  g,  g, 

0  0  0  0  0  0 

figure ( 42 ) 

subplot  (2,1,1); 

hold  on 

plot  (LINUSintensity) ; 
xxdatum=l : 1 : 512  ; 

plot (xxdatum, LINUSIntensityThreshold,  ' -r ' ) ; 
hold  off 

Title (' Intensity  vs.  Pixel  row'); 
xlabel  (' Pixel  line  number'); 
ylabel ( ' Intensity '  )  ; 

axis ( [ 0 , DATAlength, 0,2.5e5] ) ;%2e5, 8e4, 8e4, 1.5e5 
figure  ( 3 ) 

000000000000000000000 


i=l; 


%  2.2  Discard  values  outside  of  bandwidth  and  below 

threshold 

while  i  <=  DATAlength 

if  (LINUSintensity (1) <LINUSIntensityThreshold) 
LINUSintensity ( 1 ) =0 ; 

else 

LINUSintensity ( 1 ) =LINUS intensity (i)- 
LINUSIntensityThreshold; 
end 
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i=i  +  l  ; 


end 


LINUS intensity; 

DATAlength  =  length (LINUSintensity ) ; 

subplot  (3,1,3); 
plot (LINUSintensity) 

%Title  (' Intensity  vs.  Pixel  row:  Post-threshold'); 

xlabel (' Pixel  line  number'); 

ylabel ( ' Intensity '  )  ; 

axis ( [ 0 , DATAlength, 0 , 5e5 ] ) ; 


000000000000000 

0  0,0,0,0,0, 

0  0  0  0  0  0 

figure ( 42 ) 

subplot  (2,1,2); 

plot (LINUSintensity) 

Title  (' Intensity  vs.  Pixel  row:  Threshold  Subtracted'); 
xlabel (' Pixel  line  number'); 
ylabel ( ' Intensity '  )  ; 

axis ( [ 0 , DATAlength, 0,2e5] )  ;%,  ,4e4,le5 
000000000000000000000 


%  3  REINDEX  NIST 

%  3.1  Data  import  and  assignment  into  2  Nxl  matrices 

%  Graduation  of  horizontal  axes 
gradNISTx  =waveband; 

gradNISTy  =J;%%%NEW  IN  V9 

%%%  Remove  NaN  Problem  in  gradNISTs 
check  =  find(isnan(gradNISTy)); 
gradNISTy (check) =0  ; 

%Normalizing 

EilteredIntensity=EilteredIntensity/max (Eilteredintensity ) ; 
LINUSintensity=LINUSintensity/max (LINUSintensity) ; 

figure ( 4 ) 
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subplot  (2,1,1)  ; 
hold  on 

plot (NISTWavelength, Filteredintensity , 'g' ) 
hold  off 

axis ( [MinWavelength, MaxWavelength, 0,1]); 

subplot  (2,1,2) ; 
plot (LINUSintensity ) 
axis  (  [ 1 , DATAlength, 0,1]); 


000000000000000000000000000000000000000 


StretchFactor  =  Init lalStretchFactor ; 

Index  =  (FinalStretchFactor- 
Stret chFactor ) / StretchInterval+1 ; 

1  =  1; 

maxCorrMod  =  [  ]  ; 

maxCorrUnmod  =  []; 

stret chFactorAxis  =  []; 
lengthCorr  =  [  ]  ; 

figure  (5) 

while  (I  <=  Index) 

interval  =  (MaxWavelength- 
MinWavelength) / (StretchFactor*divisions)  ; 

StretchNISTx  =  MinWavelength : interval iMaxWavelength; 
StretchNISTy  = 

interpl (gradNISTx, gradNISTy, StretchNISTx) ; 

subplot (Index, 2, 1  +  2* (I-I)  ) 
hold  on 

plot (StretchNISTy, ' g ' ) %, ' . g ' ) 
plot (LINUSintensity) 

%axis ( [0,512, 0,4e5] ) ; 
axis ([0,512,0,1]); 
hold  off 

crosscorrelation  =  xcorr (LINUSintensity, StretchNISTy) ; 
maxCorrUnmod ( I )  =  max (crosscorrelation) ; 

%  For  "Pseudo-Normalization"  of  xcorr  values 

Ic  =  length (LINUSintensity) Tlength (StretchNISTy) -1 ; 

%%%  Necessary  since,  from  MATLAB  helpguide... 

%%%c  =  xcorr  (x,y)  returns  the  cross-correlation 
sequence  in  a  length  2*N-1  vector,  where  x  and  y  are 
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%%%length  N  vectors  (N>1) .  If  x  and  y  are  not  the  same 
length,  the  shorter  vector  is  zero-padded  to 
%%%the  length  of  the  longer  vector. 
lengthCorr ( I )  =  Ic; 

crosscorrelation  =  crosscorrelation  /  Ic; 
maxCorrMod ( I )  =  max (crosscorrelation) ; 
stretchFactorAxis (I)  =StretchFactor; 

subplot (Index, 2,2*1) 
plot  (crosscorrelation) 

%axis  (  [0,  1200, 0, le9] ) ; 

StretchFactor  =  StretchFactor  +  Stretchinterval ; 

1=1+1; 

end 

figure ( 6) 
subplot  (3,1,1) 

plot  ( StretchFactorAxis , maxCorrUnmod) 
subplot  (3,1,2) 

plot ( StretchFactorAxis , lengthCorr, ' r ' ) 

%  Peak  corresponds  to  point  of  maximum  correlation  growth 
%  when  referenced  to  growth  of  correlation  length 
subplot  (3,1,3) 

plot ( StretchFactorAxis , maxCorrMod) 

%  Redisplay  of  above  alone  for  emphasis 
figure  ( 7 ) 

plot ( StretchFactorAxis , maxCorrMod) 


0000000000000000000000000000 

%  Relating  pixel  line  number  to  wavelength 

[BestCorrValue, BestCorrIndice]  =  max (maxCorrMod) 

BestStretchFactor  =  Init lalStretchFactor  + (BestCorrlndice- 
1) *StretchInterval 

%BestStretchFactor  =  .168  %Manual  Override  increase 
decreases  LINUS  width 
interval  =  (MaxWavelength- 

MinWavelength) / (BestStretchFactor*divisions )  ; 

StretchNISTx  =  MinWavelength : interval : MaxWavelength; 

StretchNISTy  =  interpl (gradNISTx, gradNISTy, StretchNISTx); 
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crosscorrelat ion  =  xcorr (LINUSintensity , StretchNISTy ) ; 
maxCorrUnmodFinal  =  max (crosscorrelation); 


Ic  =  length (LINUSintensity) +length (StretchNISTy) -1 ; 
lengthCorr ( I )  =  Ic; 

crosscorrelation  =  crosscorrelation  /  Ic; 
[maxCorrModFinal , maxCorrIndice ]  =  max (crosscorrelation) 
lengthCorr  =  length (crosscorrelation) 


figure  ( 8 ) 

plot  (crosscorrelation) 

%axis  (  [0,  1200, 0, 3e9] ) ; 

lengthStretch  =  length (StretchNISTx) 

StretchNISTx ' ; 

Stretchindice  =  1 : 1 : lengthStretch; 
offset  =  maxCorrlndice- 

max ( lengthStretch, DATAlength) ; %%%FIXED ! ! ! 

%offset  =  835- 

max  (lengthStretch, DATAlength) ; %%%FIXED ! ! ! 289, sf=. 643 
pixelWidth  =  interval 

LinusMin  =  MinWavelength  -  of f set *pixelWidth 
LinusMax  =  LinusMin  +  511*pixelWidth 
LINUSwavelength  =  LinusMin : interval : LinusMax; 

figure ( 9) 

LINUSline  =  1:1:512; 

MinPlotWavelength  =  min (LINUSwavelength) 
MaxPlotWavelength  =  max (LINUSwavelength) 
plot (LINUSline, LINUSwavelength) 

axis ( [1, 512 , MinPlotWavelength, MaxPlotWavelength] ) ; 
Tit le (' Wavelength  vs.  Pixel  Row ' )  ; 
xlabel (' Pixel  Row  Number'); 
ylabel  (  'Wavelength  (Angstroms)  ' ) ; 


figure ( 10 ) 
subplot  (2,1,1); 
hold  on 

plot (NISTWave length, Elite redintensity) 
hold  off 

axis ( [MinWavelength, MaxWavelength, 0,1]); 
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Tit le (' Aligned  Pt  Standard  Data'); 
xlabel ( 'Wavelength  (Angstroms) ' ) ; 
ylabel (' Normalized  Intensity'); 

subplot  (2,1,2); 

plot (LINUSwavelength, LINUSintensity ) 
axis ( [MinWavelength, MaxWavelength, 0,1]); 

Title ( 'Aligned  LINUS  Data'); 
xlabel ( 'Wavelength  (Angstroms)  ' )  ; 
ylabel (' Normalized  Intensity'); 

%  Computation  of  spectral  properties 

SprectralRangeCamera  =  MaxPlotWavelength-MinPlotWavelength 
SpectralRangePixel  = 

SprectralRangeCamera /length (LINUS line) 

figure  (11) 
hold  on 

plot (NISTWavelength, Filteredintensity , 'g' ) 
plot (LINUSwavelength, LINUSintensity) 
hold  off 
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